home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / bisect2.f < prev    next >
Text File  |  2002-06-11  |  1KB  |  54 lines
  1. C
  2. C PAGE 77-78: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: BISECT2.FOR
  5. C
  6. C SECOND VERSION OF BISECTION METHOD (BISECT,F,G)
  7. C 
  8.       EXTERNAL F,G
  9.       DATA AF/0.0/, BF/1.0/, N/50/
  10.       DATA AG/0.5/, BG/2.0/ 
  11.       CALL BISECT(F,AF,BF,N)
  12.       CALL BISECT(G,AG,BG,N)
  13.       STOP
  14.       END 
  15.   
  16.       FUNCTION F(X) 
  17.       F = EXP(X) - 3.0*X
  18.       RETURN
  19.       END 
  20.   
  21.       FUNCTION G(T) 
  22.       G = T**3 - 2.0*SIN(T) 
  23.       RETURN
  24.       END 
  25.   
  26.       SUBROUTINE BISECT(F,A,B,N)  
  27.       FA = F(A) 
  28.       FB = F(B) 
  29.       IF( SIGN(1.0,FA) .EQ. SIGN(1.0,FB) ) THEN 
  30.         PRINT 5,A,B 
  31.         RETURN
  32.       ELSE
  33.         PRINT 3 
  34.         WIDTH = B - A 
  35.         DO 2 I = 1,N
  36.           C = A + (B - A)*0.5 
  37.           FC = F(C) 
  38.           WIDTH = WIDTH/2.0 
  39.           IF( SIGN(1.0,FA) .EQ. SIGN(1.0,FC) ) THEN 
  40.             A = C 
  41.             FA = FC 
  42.           ELSE
  43.             B = C 
  44.             FB = FC 
  45.           ENDIF   
  46.           PRINT 4,I,C,FC,WIDTH
  47.    2    CONTINUE
  48.         RETURN
  49.       ENDIF       
  50.    3  FORMAT(//3X,'STEP',10X,'C',13X,'F(C)',5X,'ERROR',//)
  51.    4  FORMAT(2X,I5,E22.14,2E10.3)     
  52.    5  FORMAT(2X,'FUNCTION HAS SAME SIGN AT',2E22.14)      
  53.       END 
  54.